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Abstract: 

This note describes a theoretical study of Wbb and Wjj final states at the Tevatron us- 
ing the NLO program MCFM. We extensively study the effect of NLO corrections with 
respect to variations of input parameters such as the minimum jet pt and the choice of 
renormalization and factorization scales. In particular, we examine possible implications 
for the Method 2 QCD background subtraction technique for it production. 
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1. Introduction 

Final states involving vector bosons accompanied by a heavy quark pair and/or light quark 
or gluon jets serve both as an arena for precision tests of QCD as well as backgrounds to 
both Standard Model and non-Standard Model physics (see for example [l]-[4]). Thus, it 
is important to utilize the most important theoretical tools available for computing their 
cross sections. In Run I at the Tevatron, predictions were obtained using a leading order 
calculation (VECBOS, [5]) supplemented by a parton shower Monte Carlo (HERWIG, [6]). 
VECBOS calculates the production of W + n jets and Z + m jets for n < 4 and m < 3, at 
leading order in the strong coupling a s . 

In Run II, there are many programs capable of calculating these processes to a much 
higher jet multiplicity than was available in Run I [7-11]. Their use at the Tevatron has 
been explored in a series of workshops at Fermilab [12]. In addition, there is now a universal 
interface (the Les Houches accord [13]) between matrix element and parton shower Monte 
Carlo programs. However, the above-mentioned matrix element programs are still leading- 
order calculations and thus have a large theoretical uncertainty. In addition, there is still 
some ambiguity in dealing with the soft and collinear cutoffs of the matrix element programs 
when matching to a parton shower Monte Carlo, and thus a risk of under or double-counting 
contributions. There are reasonable prescriptions for managing the cutoffs currently in use 
in CDF [14]; in addition, a rigorous prescription (CKKW [15]) that removes most of the 
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soft / collinear cutoff dependence in the matrix element to Monte Carlo merging has recently 
been adopted for use at the Tevatron by Steve Mrenna and Peter Richardson [16]. 

It's only at next-to-leading order (NLO), though, that the normalization of a calcula- 
tion can be taken seriously. Once more, the theoretical predictions have evolved since Run I, 
where calculations were available for W/Z+n jet production only for n < 1 (DYRAD, [17]). 
The current state of the art for calculations of this type involves W/Z plus a bb pair, or 
plus two jets [18]. MCFM provides a calculation to NLO of W/Z + bb or W/Z + jj final 
states, but at the partonic level only. Soft and collinear singularities are cancelled between 
the one-loop and tree level diagrams; b quark mass effects have still not been included in 
the NLO matrix element calculation, though, so a cut must be applied to ensure that the 
b and b are well-separated 1 . 

It is worth noting that there have been recent advances in the calculation of one-loop 
multi-partonic final states (e.g. Wbbj, Wjjj [19]) in a semi-automated fashion, but actual 
programs should not be expected for at least a year, and so will not be available for the 
first publications from Run II. It is also worth noting that the MC@NLO program [20] has 
successfully incorporated NLO matrix elements for WW, WZ, ZZ, tt and bb and Higgs in 
a parton shower Monte Carlo framework (HERWIG). Thus, the result is a fully exclusive 
final state at the hadron level. It is very time-consuming, though, to add each new process; 
luckily Wbb and Wjj remain fairly high on the priority list. 

In this note, we will examine some of the characteristics of the Wbb final states using 
MCFM at both LO and NLO. In particular, we will emphasize the impact of these calcu- 
lations on our understanding of Method 2 (which we outline in Section 5), used by CDF 
in both Run I and Run II to determine backgrounds to tt production. We will directly 
calculate the ratio Wbb/Wjj and discuss the general features of this ratio for higher jet 
multiplicities. 

2. Topologies 

The topologies for Wbb production (see Figure 1) are much simpler than those for W+2 jets 
(see Figure 2). At LO, there are basically two diagrams for Wbb production (involving only 
a qq initial state) compared to literally hundreds of diagrams for the production of W+2 jets 
(involving both qq, gq, and gg initial states). 

The bb pair invariant mass distribution peaks at low values of the mass due to the 
gluon propagator. This can be seen in Figure 3, where we show the lowest order prediction 
using MCFM. For W + 2 jet production, t-channel dijet production dominates so that the 
dijet masses will in general be higher for Wjj than for Wbb final states. Indeed, this is the 
case as can be observed in Figure 4, where the ratio of Wbb/Wjj is plotted as a function 
of the dijet mass. 

The Feynman diagram topologies may give some indication of reasonable values for 
the renormalization and factorization scales to use for the matrix element evaluation at 

1 MCFM imposes m b i > 4ml to regulate the divergence and also to ensure that the cross-section vanishes 
below the physical threshold. 
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Figure 1: Leading order diagrams for Wbb production. 
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Figure 2: A few of the leading order diagrams for Wjj production. 
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Figure 3: The bb invariant mass in Wbb events, using lowest order matrix elements. 



leading order. The next-to-leading order result may provide further evidence to support a 
particular scale(s) deemed appropriate at leading order. 
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Figure 4: The ratio of the cross section for Wbb to Wjj production as a function of the dijet mass. 
3. Partonic Cross Sections and Jet Algorithms 

MCFM is a parton-level event generator with at most 3 partons in the final state. No 
information is available at the hadron level; thus any jet algorithms must be applied to the 
1 or 2 partons that comprise any of the predicted final state jets. 

Thus far in Run II, W+jet cross sections have been measured only with cone jet 
algorithms (although the ultimate goal is to also utilize kx algorithms as well). There are 
two options: the Run I legacy algorithm (JetCLU) and the joint CDF-DO Run II algorithm 
(midpoint). The midpoint algorithm is so-named because it places a seed at the energy- 
weighted midpoint between two partons, something the JetClu algorithm does not. The 
midpoint algorithm also lacks the JetClu feature of "ratcheting" , where seed towers initially 
in the jet cone are added to the final jet energy, even if the final cone should not nominally 
include these towers. 2 Such an effect is difficult to model at the partonic level. The two 
effects end up being in the opposite direction so the differences between the two algorithms 
should be small (5% or less). One notable difference between the two algorithms is that 
the midpoint algorithm is defined in terms of the transverse momentum (pr) rather than 
the transverse energy (Er)- Details of the application of the two algorithms to partonic 
level cross sections can be found in Ref. [21]. For historical reasons, the cross sections 
generated thus far using MCFM have used a algorithm using parameters similar to 
a cone algorithm of radius 0.4; however, this should not affect any of the conclusions of 
this paper. Calculations using cone algorithms, appropriate for direct comparisons to the 
measured Run 2 cross sections, are the subject of current study [22]. 

Jet energies measured in the CDF detector have to be corrected for comparison to 
theoretical predictions. The level of correction can be tricky for comparison to calculations 

2 No seed tower left behind. 
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Figure 5: The scale dependence for LO, NLO exclusive and NLO inclusive Wbb and Wjj produc- 
tion, using our usual set of cuts. The renormalization and factorization scales are equal. 

at the NLO level. For example, one does not want to correct for the energy deposited out- 
of-cone due to perturbative gluon emission, since this is already accounted for to some level 
in the theoretical calculation. A correction should be made, however, for hadronization 
out-of-cone effects, since these are not present in the partonic calculation. The average 
hadronization correction per jet is on the order of 1 GeV. 

4. Scale Dependence 

In Figure 5, the scale dependence for Wbb and Wjj production is shown using cuts similar 
to those used in CDF [23]: 

R cone = 0A, Ai?>0.52, pj* > 15 GeV, \rj jet \<2. (4.1) 

For this calculation, the renormalization and factorization scales have been set equal. In 
perturbative QCD the freedom exists to set the two scales separately, however we have 
chosen not to do so at this stage. For inclusive cross sections, there is not the freedom 
to change the renormalization scale independently at each vertex, in contrast to a parton 
shower Monte Carlo where the scales at each vertex may be different. In Section 7, we will 
consider the impact of applying separate factorization and renormalization scales. 

In the calculations performed with MCFM, we distinguish between exclusive and in- 
clusive production, depending on whether there are exactly two, or two or more jets that 
satisfy the kinematic cuts. Contributions to three-jet final states arise only from the tree 
level Wbbj(W jjj) diagrams. 

The leading order cross sections for both processes decrease monotonically as the renor- 
malization/factorization scale increases. Both the strong coupling constant a s and the par- 
ton distribution functions (in the relevant kinematic range) decrease with increasing scale. 
At NLO, the scale dependence is reduced for both processes and for both inclusive and 
exclusive production. The logarithms that are responsible for the large variations under 
change of scale at leading order are exactly cancelled through next-to-leading order; any 
remaining scale dependence is at higher order still. 
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Figure 6: The scale dependence for LO, NLO exclusive and NLO inclusive Wbb and Wjj produc- 
tion when using a higher jet px cut of 30 GeV (upper plots) and a lower value of 7.5 GeV (below). 



At NLO, the cross section typically increases slightly as the scale decreases and then 
at some point peaks and then drops with decreasing scale, due primarily to the same 
logarithms that cancel out the scale dependence to NLO. The exact point at which the 
maximum of the cross section occurs depends both on the process under consideration as 
well as the kinematic cuts, for example the minimum jet transverse momentum. 

The scale dependence for W + 2 jets seems to be under good control for both the 
inclusive and exclusive final states, as long as a scale larger than 30 GeV is chosen. For 
exclusive Wbb production, the scale dependence is also reasonably small for scales larger 
than 30 GeV, but a considerable scale dependence remains for inclusive final states. This 
is due to the relatively large number of new channels (with gq, gg initial states) available 
for Wbbj production at NLO, and the fact that these processes are only calculated at tree 
level and thus have a leading order scale dependence. 

It is interesting to examine the impact of changing the kinematic cuts on the stability 
of the NLO calculation. In Figure 6, the study is repeated requiring the jet threshold to be 
7.5 GeV and 30 GeV. As expected, the inclusive scale dependence gets worse (better) for 
the 7.5 GeV (30 GeV) jet cut, as the 3-parton final state contribution increases (decreases). 
Note that for Wbb, the scale at which the NLO cross section becomes unstable moves to 
higher values as the jet transverse momentum cut increases. 
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Figure 7: if-factors for Wbb (solid curves) and Wjj (dotted) for exclusive final states for the 3 
different jet transverse momentum cuts. 

The K-factors (NLO/LO) are shown in Figure 7 for Wbb and Wjj exclusive final 
states for the three different jet transverse momentum cuts. The point at which the NLO 
cross section equals the LO cross section (i.e. the if-factor is 1) is relatively insensitive 
to the jet transverse momentum cuts for the case of Wjj but systematically moves out to 
higher scale values for Wbb production as the jet px cut increases. For a scale of 100 GeV 
and for a jet px cut of 15 GeV, the Wbb if-factor is 1.15 while the Wjj if-factor is 1.05. 

It is also interesting to plot the ratio of the K-factors of the two processes (Wbb, Wjj) 
as a function of the scale. The ratio of if-factors, shown in Figure 8, has a strong de- 
pendence on the jet px cut because the Wbb iT-factor does. The ratio is observed to be 
relatively constant for scales on the order of 100 GeV or above. The ratio is particularly 
unstable for low scales. 

5. Top Physics 
5.1 Two Jet Bin 

One of the most promising channels for searching for it events is a final state consisting 
of a high pt lepton plus missing transverse energy plus jets. The W boson from one of 
the tops has decayed into a an electron or muon and neutrino while the W from the other 
top has decayed into two quarks. Thus, there will be 4 partons in the final state and one 
expects it final states to have their largest contributions in the + 4 jet bins, depending 
on the jet transverse momentum cuts. Lacking NLO calculations of Wbbjj final states, 
CDF has estimated the number of non-ii events in the W+ jets sample by calculating 
the theoretical ratio of the Wbb + [n — 2) jet cross section to the W + n jet cross section 
(n = 2, 3, 4) and then multiplying this ratio by the observed number of jets in the W + n jet 
bin. Individually, such LO cross sections of course have a large scale dependence and thus a 
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Figure 8: The ratio of the if-factors for the two processes, as a function of the common renormal- 
ization and factorization scale y,. 

large uncertainty. The hope is that the .fT-factors and the scale dependence will be similar 
for the two processes and thus the LO ratio will be an accurate estimate of the NLO one. 
In CDF, this has been known as Method 2. Such assumptions can be explicitly tested in 
the Wbb/Wjj case. In Figure 8, the ratio of the K-factors for the exclusive processes Wbb 
and Wjj were plotted as a function of the scale for three different jet pt cuts. For Method 
2 to be accurate, the ratio should be near unity. This is true (K Wb ^/ Kwjj = 1-05) for a 
jet pt cut of 15 GeV, but perhaps a kinematic accident rather than a God-given truth. It 
is not true for for any scale for a pt cut of 7.5 GeV nor for small scales (20 — 30 GeV) for 
any jet px cut. 

The ratio of the Wbb/Wjj cross sections is shown in Figure 9 for both exclusive and 
inclusive production. For both cases, the ratio is examined as a function of the minimum 
cut on the jet transverse momentum. It is interesting to note that for low scales, the ratio is 
more stable at LO than at NLO. As noted earlier, this is due to the relatively large amount 
of tree-level three-parton final states that enter into the Wbb process at NLO. At NLO, the 
ratio Wbb/Wjj is approximately 1.25% for a 15 GeV jet cut, while at LO this ratio is 1.18% 
(but again very sensitive to the kinematic cuts). For lower scales, the discrepancy is much 
more extreme. At LO, the ratio has a strong dependence on the jet transverse momentum 
cut; this dependence is greatly reduced at NLO. A scale of approximately 100 GeV is in 
the region of stability at NLO. For this scale, for a jet px cut of 15 GeV, the K-factor is 
also on the order of unity (as already noted). The inclusive ratio also has the instability 
at small scales and is approximately 0.1 higher than the exclusive case for all scale values. 

The exclusive NLO ratio for Wbb/Wjj of 1.25% calculated at the scale of 100 GeV 
is larger by a factor of 1.4 than the value assumed for the 2 jet bin using ALPGEN at 
LO in the CDF Run II top analysis [23]. However, it is in good agreement with the value 
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Figure 10: The effect of a dijet mass cut on Wbb cross section and on the Wbb/Wjj ratio. 

calculated for the 2 jet bin in Run I, after an empirical factor of 1.4, determined from jet 
data [24], is applied. Thus, this analysis can be considered as a validation of the factor 
of 1.4. Note that the leading order value of the ratio in the 2 jet bin in this analysis is 
larger than the value observed in ALPGEN. This may be due in part to the b quark being 
massless in this analysis, but not in the ALPGEN calculation. 

It was noted earlier that the bb mass peaks more strongly at low mass than does the 
dijet mass. We also examine the lowest order Wbb/Wjj ratio requiring (rrijj) be 
greater than 20, 30 and 40 GeV, as shown in Figure 10. As expected, the ratio Wbb/Wjj 
decreases as the bb mass cut is increased. Furthermore, the impact of neglecting the b mass 
in the calculation decreases as the bb mass cut is increased. 

5.2 Higher Jet Multiplicity Bins and Kinematical Analyses 

NLO calculations are not yet available for the 3 and 4 jet bins; however, the inclusive ratio 
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Figure 11: The f/r distributions for Wbb(j) (lower curves) and Wjj (J) (upper curves). 

for Wbb/Wjj does give us some information about these higher jet multiplicity states. 
The topology of Wbb + (n — 2) jet final states still remains fundamentally different from 
W + n jet final states because of the importance of the g — ► bb vertex in the former. Gluon 
splitting is a negligible contribution to the W+ jets sample. The strong scale dependence 
observed for Wbb + 2 jets at NLO may lead to some wariness about predictions for the 
behavior of higher jet multiplicity states. However, as we have argued, the scale sensitivity 
results from the relative importance of the Wbb +jet tree level diagrams; this will not 
happen at NLO for the 3 or 4 jet bin. In lieu of NLO higher jet multiplicity calculations, 
we can try to form observables from suitably inclusive quantities for which the NLO Wbb 
and Wjj calculations may form an adequate approximation. One such observable is Ht] 
here, Ht is defined as the sum of the transverse momenta of all of the jets, the lepton and 
the missing pt in the event. The Ht distributions for Wbb(J) and Wjj(j) are plotted in 
Figure 11. The (LO) distributions for Wbbj and Wjjj naturally peak at at a higher value 
of Ht due to the requirement of an additional jet. 

For Wjj the Ht distributions are similar at LO and NLO. The Wbb distribution 
is steeper at LO than at NLO. In Figure 12, the cross sections for Wbb and Wjj are 
normalized for shape comparison. There it can be clearly observed that, although, the 
Wbb and Wjj shapes are similar at LO, the Wjj distribution is steeper than that of Wbb 
at NLO. Any assumption that the shapes of the two processes are similar (as for example 
in some kinematical analyses) can be dangerous. 

In Figure 13, the ratio of Wbb to Wjj, at LO and NLO (both inclusive and exclusive), 
as well as the (LO) ratio of Wbbj to Wjjj is plotted as a function of Ht- Here, the 
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Figure 12: The Ht distributions for Wbb(j) and Wjj(j), normalized to the same area. 

similarity of the Wbb and Wjj distributions at LO can be clearly observed, as well as the 
fact that the NLO ratio of Wbb to Wjj increases by more than a factor of two over the 
Ht range of the plot (which is approximately the same as the Ht range used in the top 
analysis). Similar results to those above have also been obtained when replacing Ht by 
the scalar sum of the jet transverse energies. 

Given the shape differences between LO and NLO observed for global variables such 
as these, it is important to discover if the effects are the same for less inclusive observables. 
In Fig. 14 the ratio of Wbb to Wjj is plotted as a function of the lead jet pt- In contrast 
to the previous cases, the ratio is flat at NLO for both inclusive and exclusive final states. 
Thus, this is an example of a variable which appears to be safe with regard to a change of 
shape at NLO. 

In Figs. 15 and 16 we show the relative contributions to the NLO Ht and largest jet 
Pt cross sections of the Wbb(Wjj) and Wbbj(Wjjj) subprocesses. In the NLO inclusive 
results, the contribution to the Ht distribution for Wbbj and Wjjj events is negligible at 
small Ht and dominant at large Ht, leading to the change in shape seen in Fig. 13. This 
can be contrasted with the extra jet contribution to the largest jet pt distribution, which 
is never dominant over the range shown. 



6. Two Scale Studies 



Thus far, we have set the renormalization and factorization scales to the same value. A 
version of MCFM has been modified to allow the separation of the two scales. In Figure 17, 
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Figure 13: The ratio of the Ht distributions for Wbb(j) and Wjj(j). 

the Wjj cross section has been plotted with the two scales set the same (both at LO and 
NLO) as well as for one scale pre-factor being set to the inverse of the other (at NLO). A 
similar plot is shown in Figure 18 for the case of Wbb. As expected, the scale dependence 
is large at leading order; at NLO the cross section is reasonably stable for all of the scale 
choices plotted. Thus, the ratio of the two cross sections used in Method 2 should also be 
stable at NLO for the case of unequal renormalization and factorization scales and there 
is no anomalous enhancement of the Wbb cross section at low renormalization scales. 

7. Other Theoretical Uncertainties 

One measure of possible NNLO corrections is supplied by the scale dependence that remains 
at NLO. This was shown in Figures 4-9 where it was observed that the scale sensitivity 
is under control for choice of scale larger than 50 GeV or so and for reasonably large jet 
Et cuts. The ratio can also be sensitive to the choice of parton distribution functions 
(PDF's). We have used the 40 error PDF's supplied with the CTEQ6M central fit PDF to 
estimate the contribution of PDF uncertainty to the ratio of Wbb to Wjj production. As 
noted earlier, different initial states contribute to the two different processes, so that PDF 
uncertainties will not necessarily cancel. This is the first time (that we know of) that such 
a study has been carried out for Wbb/Wjj final states. In order to more easily calculate 
the PDF uncertainties, we have used the LHAPDF interface provided in the most recent 
version of MCFM [25] . With this version, the parton luminosities for each PDF member 
can be calculated very quickly at each Monte Carlo integration point. Thus, the cross 
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Figure 14: The prjeti distributions for Wbb(j) and Wjj(j). 

sections using all members of a given PDF set can be calculated in one Monte Carlo run. 
The integration is weighted by the central PDF luminosity and the result for each error 
set is recorded separately. 

The PDF uncertainties thus calculated are shown in Figs. 19 and 20 for the Ht and 
largest jet px variables. The PDF uncertainties are reasonably small. 

8. Conclusions and Future Plans 

Method 2 has proven to be a useful tool for CDF analyses both in Run 1 and in Run 
2. In this paper, we have shown that Method 2 holds its validity at NLO for non- 
inclusive variables, for jet Et cuts of the order of 15 GeV or greater, and for renormaliza- 
tion/factorization scales of the order of the W mass. 

In an ideal world, we would have available NLO calculations of Wjjj, Wjjjj, Wbbj 
and Wbbjj. Since their availability within the next year is unlikely, we will have to continue 
to rely upon LO predictions of these final states. In a future paper [22], we will attempt to 
use the NLO processes within MCFM to directly compare MCFM predictions to (1) Wjj 
and Wbb observables in CDF Run 2 data and (2) to enhanced LO predictions (for example 
using the CKKW scheme). Good agreement with the latter for Wjj and Wbb final states 
may give some confidence in the extrapolation to higher jet multiplicity final states. 3 In 
particular, it is believed that some explicit higher order corrections are incorporated into 

3 We will also consider some of the subtleties generated by comparing parton level calculations to data 
cross sections measured with a relatively small (0.4) cone size. 
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Figure 15: The relative contributions to the NLO Ht cross section of the Wbb(Wjj) and 
Wbbj(Wjjj) subprocesses. 

the Sudakov form factors through the use of the CKKW procedure. We can try to explicitly 
test this. 
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Figure 17: The scale dependence of the Wjj exclusive cross-section, using our usual choice of 
varying renormalization and factorization scales together (fip — [ir — RMyy, for 1/2 < R < 2), as 
well as the choice of varying them in opposite directions (fxp — RMyy, = 1/RMyy). 
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Figure 18: The scale dependence of the Wbb exclusive cross-section, using our usual choice of 
varying renormalization and factorization scales together (fip = fin, = RMw, for 1/2 < i? < 2), as 
well as the choice of varying them in opposite directions (/if — RMw, = 1/RMyy)- 
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Figure 19: The PDF uncertainty for the ratio of Wbb and Wjj, plotted as a function of Ht, 
calculated using the CTEQ6 error PDF set. 
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Figure 20: The PDF uncertainty for the ratio of Wbb and Wjj, plotted as a function of the sum 
of the jet transverse energies, calculated using the CTEQ6 error PDF set. 
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